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Abstract 

For the example of the infinitely deep well potential, we point out some paradoxes which 
are solved by a careful analysis of what is a truly self- adjoint operator. We then describe 
the self-adjoint extensions and their spectra for the momentum and the Hamiltonian 
operators in different physical situations. Some consequences are worked out, which could 
lead to experimental checks. 
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1 Introduction 

In most French universities, quantum mechanics is usually taught in the third year courses, 
separately from its applications to atomic, molecular and subnuclear physics, which are dealt 
with during the fourth year. In such "first contact" lectures, many mathematical subtleties 
are necessarily left aside. However, even in such commonly used examples as infinitely deep 
potential wells, overlooking the mathematical problems leads to contradictions which may be 
detected by a careful student and which have to do with a precise definition of the "observables" 
i.e. the self-adjoint operators. 

Of course, experts in the mathematical theory of unbounded operators in Hilbert spaces 
know the correct answer to these questions, but we think it could be useful to popularize these 
concepts among the teaching community and the more mature students of fourth year courses. 
In particular, the role of the boundary conditions that lead to self-adjoint operators is missed 
in most of the available textbooks, the one by Ballentine || p. 11] being a notable exception 
as it includes a discussion of the momentum operator. But there, we find only two references 
relevant to the subject. The first one considers a particular self-adjoint extension of the 
momentum and of the Hamiltonian for a particle in a box, which is interpreted as describing 
a situation with spontaneous symmetry breaking. The second one M mentions the self-adjoint 
extensions of the Hamiltonian for a particle in a semi-axis and its relevance, first pointed out 
by Jackiw ||, to the renormalization of the two dimensional delta potential. 

The aim of this paper is to emphasize the importance of the boundary conditions in the 
proper definition of an operator and to make available to an audience of physicists basic results 
which are not so easily extracted from the large amount of mathematical literature on the 
subject. 

The paper is organised as follows : in Section 2 we discuss some paradoxes met in the study 
of the infinite potential well. Then, in Section 3, we present a first analysis of the boundary 
conditions for the self-adjoint extensions of the momentum operator. 

In Section 4 we introduce the concept of deficiency indices and state von Neumann's theorem. 
In Section 5 we apply it to the self-adjoint extensions of the momentum operator for which the 
spectra, the eigenf unctions and some physical consequences of these are given. We hope that, 
despite some technicalities needed for precision (which can be omitted in a first reading), the 
results are of easy access. The reader interested in these technical aspects may consult the 



references [0] and [fL2 



Then, in section 6, we describe the self-adjoint extensions of the Hamiltonian operator 
in various settings (on the real axis, on the positive semi-axis and in a box). Several physical 
implications are analysed, while in section 7, we use different constraints from physics to reduce 
the set of all possible self-adjoint extensions. 

We have gathered in appendix A some technical details on the extensions of the momentum 
operator and in appendix B we discuss the spectra of the Hamiltonian operator for a particle 
in a box. A proof for parity preserving self-adjoint extension is given in Appendix C. 

2 The infinite potential well : paradoxes 

Let us consider the standard problem (see for example |L0], p. 299] or || p. 109]) of a particule 
of mass m in a one dimensional, infinitely deep, potential well of width L : 

V(x) = 0, xG]-J,+i; V(x)=oo, \x\>^. (1) 
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Stationary states are obtained through the Schrodinger (eigenvalue) equation 

H(j){x) = E(j){x) 

and the vanishing of their wave function at both ends. This means that the action of the ha- 
miltonian operator for a free particle, unbounded on the closed interval [— +■§], is defined 
by: 

H = -^D\ V(H)= e £ 2 (-§,+§) , 0(±§) = o}, (2) 

where D is the differential operator — and T>(H) is the definition domain of the operator H. 

dx 

Two series of normalised eigenfunctions of opposite parity are obtained. They vanish outside 
the well and for x E [—-§,+§] they write : 

[2 . 2n7ix ti 2 f2nix^ 2 

odd ones : $n{x) = \ ~F sm — F — I > = o — 



L L " 2m V L 

T . . /2~ .(271-1)7^, „, /I 2 /(2n-l)7rV ^ 
ewer* ones : tf n (a;) = y^cosP ^— ] , E' n = — l± . (3) 

where n is a strictly positive integer. The functions $ n (x) and are continuous at x — ±-| 

where they vanish. 

A question of fundamental importance arises : is the Hamiltonian operator H a truly self- 
adjoint operator ? To discuss more thoroughly this question let us consider a particle in the 
state defined by the even, normalised wave function : 

= -yl^ 2 - x } ' |x| - \ ; = ' |x| - \ ■ (4) 

It may be expanded [17] on the complete basis of eigen functions of H given in (|3]) : 

u n v n \x) , u n = yw n , — - 

n=l 

Let us define also, for further use 



¥(*) = 5>*»(*) , b n = (* n> v]/) = ^ ^ . (5) 



= ~—D 2 m{x) = -J^r , -L/2 < x < +L/2, (6) 

and let us begin with some elementary computations : the mean value of the energy and its 
mean-square deviation in the state (f|). On the one hand we have 

<E>-}^ \b n \ E n - ^ _ - mL2 , (7) 

n=l n=l v ' 

but on the other hand 

x , ~x 30£ 2 f +L ' 2 , 2 L 2 1 5h 2 10 , 

< E >= (*, OT ) = (*, *) = [x 2 - -]dx = ^ = - 2 E[. 
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These results are coherent. Things are different for the energy mean-square fluctuation. On 
the one hand 



240fr 4 A 1 30fr 4 

WW = - 

leads to 



m %2£4 ^ f 2n - 1)2 m 2 L 4^ 
n=l v ' 



, fc2 

AE = V< E 2 > — < E > 2 = V5—^, (9) 

and on the other hand 

< E 2 >= # 2 ^) = #\I>) = !! 

In order to understand the origin of the paradox, let us come back to the definitions. The 
probability of being in the eigenstate <p n of energy e n being given by \l/)| 2 , one obtains 

oo oo oo 

< E 2 >= J24\(<Pn,n 2 = £)<£(0n,*)OMn) = J2Wn,^,H<p n ) 

n=l n=l n=l 

where the reality of the eigenvalues of the Hamiltonian has been used. If H were self-adjoint, 
one would obtain with the help of the closedness relation 

°° ~ 30h 4 

< E 2 >= y70„, H*)(HV, cj> n ) = (HV, HV) = (*, *) = —— . (10) 

n=l 



in agreement with the direct calculation (|9|). But, if the self-adjointness of H was used once 
more, one would get 

< E 2 >= (HV, HV) = (*, H 2 ^) = (11) 

which is necessarily wrong. In fact, in QIDD, we used (correctly, as shown by the standard proof 
using an integration by parts ) the self-adjointness of H when it acts in the set of functions 
that vanish at both end-points of the well 



(H<p n , tt) = (0 n , HV) , H<f) n ) = (HV, (j> n ) ; 

on the contrary, in ([0]), the function \I> does not belong to that set and, consequently, in the 
integration by parts, the integrated term remains and 



These simple calculations show that the problem lies in the definition of the action of the 
operator if on a function \I/ that does not vanish at the end-points. 

To summarise, we came up against the difficulty of the definition of a self-adjoint operator 
in a closed interval [-L/2, +L/2] as an extension of a differential operator —-^D 2 , question 
already solved by mathematicians in the thirties. Before explaining this theory in a simple 
manner, we analyse in the next Section the momentum operator —ihD. 
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3 Self-adjoint extensions of the momentum operator : 
a first approach 

Let us consider the one-dimensional momentum operator P = —ihD in a closed x interval. 
Let us take for domain T> the following space 

V(P) = {0, 0' G £ 2 ([0, L\) ; 0(0) = 0(L) = 0} . 

The vanishing of 



(-0, —ihDcj)) — {—ihDip, 0) = I dx 



L 



i){x){-ih—^-) - (ih—j—)(j){x) 

= -ih I dx-£-${x)<f>(x)] = -ih$(L)4>(L) -^(0)0(0)] . (12) 

implies that P is a symmetric operator in T> . But P is not a self adjoint operator even if its ad- 
joint pt = —ihD has the same formal expression, but it acts on a different space of functions. 
Indeed, 

P(P') = {-0, -0' £ £ 2 ([0, L\) ; no other restriction on ip(x)} . 

With (|12|) , one easily sees that the adjoint of the operator P\ = —ihD acting on the subspace 
of £ 2 ([0, L]) such as 

0(L) = A0(O) , where A G C 

is the operator Py where A' = 1/A. As a consequence, a candidate family of self-adjoint 
extensions of the operator —ihD, depending on a complex parameter A = 1/A , i.e. a phase 
A = e id , 9 G [0, 2tt] is: 

P e 0(x) = -ihD<f>{x) , V(P g ) = {0, 0' G £ 2 ([0, L}), 0(L) = e' e 0(O) } (13) 

Notice that for 8 = 0, one recovers the usual periodic boundary conditions. 

Conclusion : A symmetric differential operator acting on a given functional space is not 
automatically a self-adjoint operator and may have none, a unique or an infinity of self-adjoint 
extensions. In the next Section, we give some mathematical results on the theory of self-adjoint 
extensions of a differential operator in a Hilbert space and "deficiency indices" . 



4 Deficiency indices and von Neumann's theorem 

Since this Section makes use of mathematical terminology, let us begin with some precise 
definitions. 

Let us consider a Hilbert space H. An operator (A,V(A)) defined on 7i is said to be 
densely defined if the subset T>(A) is dense in 7i, i.e. that for any -0 G Ti one can find in 
T>(A) a sequence 0„ which converges in norm to ip. 

An operator (A,T>(A)) is said to be closed if n is a sequence in T>(A) such that 

lim n = 0, lim A<p n = ip, 

then G V(A) and A(fi = ip. 



6 



Bonneau, Faraut, Valent 



Let us recall the definition of the adjoint operator of a (in general not bounded) operator 
H with dense domain T>(H). The domain T>(W) is the space of functions ip such that the 
linear form 

<t> — > 

is continuous for the norm of Ti. Hence there exists a ^ G Ti such that 

(^,i/0) = (^,0). 

One defines Wtp = i]v . A useful result is that the adjoint of any densely defined operator is 
closed, see [|I|, p. 80, vol. 1]. 

An operator (H, T>(H)) is said to be symmetric if for all 0, ip G T>(H) we have 

(H</>,il>) = (</>, Hif>). 

If T>{H) is dense, it amounts to saying that (H\V(H)) is an extension of (H,T>(H)). 

The operator with dense domain T>(H) is said to be self-adjoint if D(W) = T>(H) and 
#t = if. 

In this Section we will assume that (A, T>(A)) is densely defined, symmetric and closed and 
let (A\V(A*)) be its adjoint. 

One defines the deficiency subspaces Af± by 

N+ = {if> e V(Ai) , Aty = z + i/>, Imz + >0}, 
JV_ = G P(At) , Aty = Z-t/), Imz_<0}, 

with respective dimensions n + ,n_. These are called the deficiency indices of the operator A 
and will be denoted by the ordered pair 

The crucial point is that n + (resp. n_) is completely independent of the choice of z + (resp. 
z_) as far as it lies in the upper (resp. lower) half-plane. It follows that a simple way to 
determine (n + , n_) is to take z + = iX and z_ = —iX with an arbitrary strictly positive constant 
A needed for dimensional reasons. 



The following theorem, first discovered by Weyl [16] in 1910 for second order differential 
operators and generalized by von Neumann |15| in 1929, is of primary importance 

Theorem 1 For an operator A with deficiency indices (ra + ,n_) there are three possibilities : 

1. If n+ = n_ = 0, then A is self-adjoint (in fact this is a necessary and sufficient 
condition). 

2. Ifn + — n_ = n > 1, then A has infinitely many self-adjoint extensions, parametrized 
by a unitary n x n matrix (i. e. n 2 real parameters) . 

3. If n + 7^ n_, then A has no self-adjoint extension. 

The application of this theorem to differential operators requires still a lot of work : even if 
we start from an operator P which is formally self-adjoint, this does not prove that P is truly 
self-adjoint because the domains T>(P) and D(P<) will be different in general. 

For a given differential operator P one has to solve three problems : 

1. Find a domain T>(P) for which the formally self-adjoint operator P is symmetric and 
closed. 
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2. Compute its adjoint (P* ,T)(P')) and determine the deficiency indices of PL 

3. When they do exist, describe the domains of all the self-adjoint extensions. 

A whole body of theory has been built up to solve these problems and is given in many 
text-books (for instance P|,[^2| )■ In the next Section we describe the results for the simplest 



case of the momentum operator P = —ihD, referring for the proofs to |T|, vol. 1, p. 106-111]. 



5 Self-adjoint extensions of the momentum operator 

Let us apply the previous analysis to the momentum operator P = —ihD , in three different 
"physical" situations : first on the whole real axis and in this case we conclude to a unique 
self-adjoint extension, second on the positive semi-axis and in this case there is no self-adjoint 
extension, and third in a finite interval [0, L] in which case there are infinitely many self-adjoint 
extensions, parametrized by U(l) , i.e. a phase. The momentum operator is certainly the 
simplest differential operator to begin with and it already exhibits all the possibilities described 
in von Neumann's theorem. For each physical situation corresponding to position space being 
some interval (a,b), finite or not, the maximal domain on which the operator P = —ihD has 
a well defined action will be called T> max (a, b). In this Section, we apply the previous theorem, 
postponing some mathematical details to the Appendix A. 

Let us consider the Hilbert space TC = L 2 (a, b) and to use von Neumann's theorem, we have 
to determine the functions ip±(x) given by 

h 

P^tp±(x) = —ihDijj±(x) = ±i—Ajj±(x). 

a 

For dimensional reasons we have introduced the constant d > 0, homogeneous to some length. 

An easy integration gives ip±(x) = C±e Tx ^ d . Then we have to discuss the different intervals 
(a, b). 



5.1 The operator P on the whole real axis 

None of the functions ip±(x) belong to the Hilbert space L 2 (M) and therefore the deficiency 
indices are (0,0). Hence we conclude that the operator (P, D max (IR)) is indeed self-adjoint, 
in agreement with the heuristic considerations given in the standard textbooks on quantum 
mechanics. Moreover, the spectrum of P on the real axis is continuous, with no eigenvalues. 



5.2 The operator P on the positive semi-axis 

Among the functions ip±(x) , only ip + belongs to L 2 (0, +oo) . We conclude to the deficiency 
indices (1,0) and therefore, by the von Neumann theorem, P has no self-adjoint extension. 
This is a fairly surprising conclusion, since it implies that the momentum is not a measurable 
quantity in that situation ! 



5.3 The operator P on a finite interval 

Since we are working on a finite interval, both ip±( x ) = C±e Tx ^ d belong to L 2 (0,L) and the 
deficiency indices are (1,1). 
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From von Neumann's theorem, we know that the self-adjoint extensions are parametrized 
by U(l) , i.e. a phase e l9 , in agreement with the result of section 3. Denoting these extensions 
by Pq = (P, Vg), they are given by 

T> e = {i>eV mBX (0,L), t/)(L) = e i6 i/)(0)}, 9e[0,2ir]. (14) 



Moreover, the spectra are purely discrete. Using the boundary condition ( |14]) , the eigenval- 
ues and eigenf unctions are easily shown to be 

P <j) n (x,9) = ^z/0 n (x,6O, v = n + ^-, n = 0, ±1, ±2, . . . 

1 r_. xi ( 15 ) 

(0m; 0n) 3mn- 



IX, (7 I = — ^= cxp 

v 7 ! 



2mv— 



As the phase 9 appears in the eigenfunctions any measurement of the momentum of a given 
system should, in general, depend on it. To display this, let us go back to the state (|j). After 
a translation, we are left with the wave function 

, W .^(L-,). 

Its eigenfunction expansion is 

n=+oo 
n=— oo 

with coefficients 

. /on r r,;^ ^ q /o^ i 

e" i9/2 , for ^ 0, (16) 



and 



. sin(0/2) 
cos(#/2) - 



7TU 

c = — - , c n = - n = l,2, ... for 9 = 0. (17) 

6 zir^n 2 

So the probability to find the particle with a momentum — - — , being equal to |c n (6 l )| 2 , is really 

9 dependent. Of course one would like to have a physical argument which gives some preferred 
value of 9. 

Let us conclude with the following remarks : 

1. The textbooks which do study the momentum operator in a box (0 and |], vol. 2, p. 
1202]), usually consider (using physical arguments) only the self-adjoint extension correspond- 
ing to the periodic boundary condition (i.e. 9 = 0) which is certainly the simplest (but still 
arbitrary) choice. The anti-periodic boundary condition (i.e. 9 = n) has been considered by 
Capri in ||. 

2. For a particle in a box, it is often argued that the "physical" wave function should 
continuously vanish on the walls x = and x = L, ensuring that the presence probability 
vanishes continuously for x < and for x > L. One should realize that the continuity of the 
measurable quantity 

Pr(0<x<w)=/ \(p(x)\ 2 dx, ue[0,L] 

1 1\ 
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is ensured as soon as the integral / \4>{x)\ 2 dx does converge and does not require any con- 

Jo 

tinuity property of 4>(x). Specializing this remark to the eigenfunctions of Pg we observe that 
\4> n (x,8)\ 2 does not vanish continuously at x = but nevertheless the physical quantity 

Pr(0 < x < u) = j 

vanishes continuously as it should, for u — > 0. 

3. The existence of normalisable eigenfunctions of the momentum operator has an impor- 
tant consequence : the Heisenberg inequality AX ■ AP > h/2 no longer holds. Indeed, for the 
state 4> n (x,9) given in relation ([i5|) , one has AP = and AX = L/2. On the contrary, on 
the whole real axis the spectrum is fully continuous (no normalisable eigenfunctions), and the 
momentum probabilities are related to the Fourier transformed wave function. As the widths 
in x-space and in p-space are inversely proportional, the Heisenberg inequality follows. 

4. If one identifies the variable x with the angular variable p G [0, 2n] of polar coor- 
dinates, then the angular momentum is L z = —ih—. The previous remark shows that the 

dip 

inequality Aip ■ AL Z > h/2 can be violated, even by wave functions periodic in the angle ip. 



6 Self-adjoint extensions of the Hamiltonian 

In the same setting as in the previous section, we consider now the Hamiltonian operator 
H = —D 2 . We work in the Hilbert space L 2 (a, b). The maximal domain in which the operator 
D 2 is defined will again be called D max (a, b). To compute the deficiency indices we solve 

-D 2 <p(x) = ±ik%<f>{x), k >0, (18) 

and get 

<j) ± = a±e k±x + b±e~ k±x , fc± = ^^-k . (19) 

a/2 



6.1 The Hamiltonian on the whole real axis 

The physical situation corresponds to a free particle moving in a one dimensional space. The 
Hilbert space is 7i = L 2 (1R) which implies <fi± ^ 7i and the deficiency indices (0,0). It follows 
that on the real axis there is a unique self-adjoint extension of the Hamiltonian, with a fully 
continuous spectrum, in full agreement with the physicist understanding of this case. 



6.2 The Hamiltonian on the positive semi-axis 

The physical problem is that of a free particle in front of an infinitely high wall for x < 0. In 
the Hilbert space Ji = L 2 (0, +oo) we have the solutions to equation (|18|) given by 

4>± = b±e- kox/V2 e ±ikox/V2 , 

leading to the deficiency indices (1, 1), and therefore to infinitely many self-adjoint extensions 
parametrized by £7(1). 
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The corresponding boundary conditions are 

(0'(O) - £0(0)) = e iC V(0) + £0(0)), a e [0, 2tt], 

which are equivalent to 

0(O) = A0'(O), A = -tan(a/2), AelU{oo}, (20) 

see []I], vol. 2, p. 187, 204]. The boundary condition 0'(O) = corresponds to A = oo . Physicists 
use the particular extension with A = 0, see for instance [|K], p. 328] and |14], p. 33]. 



Let us now discuss the energy-spectra of a particle confined in the region x > . When 
the particle energy E is positive, we can compute the reflexion coefficient for this infinitely 
high barrier in order to compare the predictions given by the different extensions. The wave 
function is 

h 2 k 2 

Mx) = Ae~ ikx + Be ikx } E = , k > 0. (21) 

2m 

Let us define the reflection amplitude and reflection probability by 

r(*) = ~ R(k) = \r(k)\ 2 . 

Imposing the boundary condition (|2"0"|) we get 



1 + iXk , , 

r(k) = — — =► R = l. 22 

1 — i\k 

Remarkably enough the physical content (i.e. R = 1 !) of all the extensions is the same : the 
wall acts as a perfect reflector. 

This is not quite true for the bound states 

E = ~^> P>°> <Kx)=Ae-", 

for which ( pUl) implies (1 + \p)A = 0. There will be a bound state with p = — 1/A only for 
A < and different from oo. Its energy and normalised wave function are 

E = ~^ A <°> <K*) = \\^*-* m - ( 23 ) 

As far as an infinitely high wall is feasible experimentally, the existence (or non-existence) of 
this negative energy will act as a selector of some self-adjoint extensions. 

If experiment rules out the negative energy state, or if one is reluctant to accept negative 
energies for the Hamiltonian, there are still many possible extensions, with A > or A = oo. 

In an attempt to lift this degeneracy, we consider the simplified deuteron theory described 
by the potential 

oo for x < 0, 
V(x) = { -V for < x < a, V > 0, (24) 
for x > a. 
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The wave function is well known to be 

x < a : <fii(x) = A sin kx + B cos kx, 
x > a : 02 0*0 = C e~ px , 



E + V 
E 



h 2 k 2 

2m ' 
2m 



fc > 0, 
p > 0. 



(25) 



We next impose the boundary condition ( p0|) and the usual continuity conditions at x = a. 
Using the notations X = ka and K = pa , we get that the bound state energy is given by the 
solution of the system 



A > 



Y = -X- 



L - (A/a)X tan X 
tan X + (X/ajX ' 



LE| = + ( y 



(26) 



In the case of the deuteron, the absolute value of the binding energy \E\, is roughly equal to 
2.2 MeV. Its size is a = 2 F and we take 2m = M where M is the nucleon mass. It follows 
that Y = 0.46. For a given value of A, we have to solve for X, and then recover the potential 
Vq. Numerical analysis gives the following dependence on Vq with respect to the parameter A : 



A/a 





0.1 


0.2 


0.5 


1 


2 


5 


10 


100 


oo 


V (MeV) 


36.8 


31.5 


27.5 


20.5 


15.3 


11.5 


8.59 


7.50 


6.47 


6.34 



Let us observe that the parameter A, describing the different extensions, does indeed have 
an effect on physical quantities (as already observed for the momentum operator, in subsection 
5.3) and in fact experiment, not just theoretical prejudices, should decide which is the "right" 
value for it. 



6.3 The Hamiltonian on a finite interval 

This last case corresponds to a particle in a box : 16 [0, L\. From a mathematical standpoint 
the situation is quite similar to the one already experienced with the momentum operator in 
the previous section, but up to our knowledge, it did not appear before in the literature. So we 
give some details in the main text. 

One starts from the operator (H,Vq(H)) such that 

V {H) = {</>€ P max (0, L) and 0(0) = 0(L) = 0'(O) = 0'(L) = 0.} 

It is densely defined and closed, with adjoint 

H* = H, V(H^)=V max (0 } L). 

Since all the solutions of equation (|i~8"D belong to L 2 (0, L), the deficiency indices are now (2, 2) 
and the self-adjoint extensions are parametrized by a £7(2) matrix. 

To describe these self-adjoint extensions, it is natural to introduce the sesquilinear form, for 
4> and ip in £> max (0,L), 

S(0,^) = ^((^V,^)-(0,^V)) (27) 
which depends only on the boundary values of and ip. Specializing to ip = we have 

B(<P, 0) = 1 U'(L)W) ~ <P(L)Wj ~ 0'(O)0(O) + 0(O)0 7 (O)) . (28) 
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The identity 

(xy -xy) = ^(\x + iy\ 2 - \x - iy\ 2 ) , (29) 

applied to x — L<j)'(L), y = <j)(L) and x = L0'(O), y = 0(0) brings relation (|28|) to 

4LB(4>, 0) = \L<f>'(0) - £0(0) | 2 + \L(f>'{L) + i<P{L) | 2 - |L0'(O) + #(0) | 2 - \L<p'(L) - i<f>(L) | 2 . (30) 

The domain of a self-adjoint extension is a maximal subspace of "D max (0, 1/) on which the 
form B(<f),<j)) vanishes identically. These self-adjoint extensions are parametrized by a unitary 
matrix U, and will be denoted Hjj = (H,T>(U)), in which T>(U) is the space of functions 
in T> max (0, L) satisfying the following boundary conditions 

L0'(O) - £0(0) \ =u f W IP) + »W) \ 
L(p'(L) +i(p(L) J \L<f>'{L)-i<f>{L) J' 

Notice the arbitrariness in the choice of the ordering of the coordinates L(f>'(0) ± £0(0) and 
L(j)'(L) =F i(f)(L). The crucial observation is that whatever the choice of coordinates is, the 
arbitrariness of the self-adjoint extensions remains described by a U (2) matrix. 

These boundary conditions describe all the self-adjoint extensions H v = (H,V(U)) of a 
particle in a box. Moreover, thanks to the useful theorem, proved in [|1^, vol. 2, p. 90], stating 
that for a differential operator of order n with deficiency indices (n, n) all of its self-adjoint 
extensions have a discrete spectrum, we know that all the spectra of the Hjj are fully discrete. 
Leaving the details of these spectra to the Appendix B, we only give the results. 

Let us parametrize the unitary matrix U as : 

U = e** M, det M = 1, => det U = e 2 ^, $ e [0, vr] (32) 

where M is an element of SU(2), i.e. a unitary matrix of determinant 1. The range of i/j is 
restricted to n instead of 2n because the couples (0, M) and (0 + 7r, —M) give rise to the 
same unitary matrix U. Notice also that it follows that the points ip = and if) = it are to 
be identified. 

To parametrize the matrix M , we used the Pauli matrices 






and the notation : n ■ f = n\ T\ + ri2 Ti + t% . With coordinates m = (mo, m) constrained 
by 

ml + m ■ m = 1 <^=^> m G § 3 , (33) 



M writes 



mo — «m 3 — m 2 — imi 
M= | " ) =rn I -im-r. (34) 

m 2 — imi mo + im.3 



Then, starting from the boundary conditions (|31~D, we obtain the spectra for the Hamiltonian 
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in a box (see details in Appendix B) : 

a) E = ^ > : 

2s[sin'0coss — mi] = sin s [cos ip(s 2 + 1) — m (s 2 — 1)], 

b) E = : s -> in result (35 - a) <^> 

2 sin ^ — cos ip = 2mi + mo, 

c) E = — — < : s = zr in result ([35]— a) <^> 

L 2 



2r[sinV coshr — mi] = sinhr[— cos ip(r 2 — 1) + mo(r 2 + 1)]. 



Remarks 



(35) 



1. The eigenvalue equations are independent of the parameters (rri2, ms). As shown in 
appendix B, this follows from their invariance under the transformation 

M -> M' = e- dTl/2i Me +dTl/2i . 



Let us point out that this invariance is specific of the spectra, not of the eigenfunctions. 

2. The existence of negative energies seems rather surprising since P 2 = —D 2 is a formally 
positive operator. That this is not generally true can be seen by computing 

(0, H4>) - (P(j>, P(j>) = 0(O)0'(O) - 0(L)0'(L), G D v . 

If the right hand side of this relation is positive, then the spectrum will be positive, an issue 
which depends on the extension Hu considered (see section 7.3). 



7 Restrictions from physics on the self-adjoint exten- 
sions 

In the previous section we have described all the possible self-adjoint extensions of the operator 
Hu as they follow from operator theory. Now we examine which extensions are likely to play 
an interesting role according to arguments from physics. 



7.1 Extensions preserving time reversal 

Let *$?(x, t) be a solution of the Schrodinger equation 

= -^^(M) (36) 

inside the box. The time reversal invariance of this equation means that if ty(x,t) is a solution 
of (|36"D , then *&(x, t) is also a solution. If we consider a stationary state of definite energy E 
with the wave function 

#(z,t) = (j) E (x)e- i ^, 
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the previous statement implies that 4>e{x) and <Pe{ x ) are two eigenfunctions of the Hamiltonian 
H with the same eigenvalue E. One can therefore choose real eigenfunctions by taking the linear 
combination 4>e{. x ) + 4>e( x )- 

The shortcoming in this argument is that the boundary conditions (|3l|) do not lead neces- 
sarily to real eigenfunctions 4>e{ x )- Among all of the self-adjoint extensions of the Hamiltonian 
only some subclass will have real eigenfunctions. These extensions will be said to be time 
reversal invariant. 

To determine all of these extensions, we merely observe that, using the notations 

i>±{%) — L<p'(x) ± i<fi(x), 
the reality of 4>(x) implies ip±(x) = tp T (x). Taking the complex conjugate of relation (|3T| ) gives 

M0)\ V (*-V)\ v V (*+V>\ (37) 

) \ ML) J \ V-(L) J 

Since ip+(0) an d ^-(-^) cannot vanish simultaneously, we conclude to 

det(l -UU) = 0. (38) 



Using for U the coordinates given by (p4|), easy computations give m 2 = and, correspond- 
ingly, the matrix 

mo — —imi \ 

U = e # ( with V e [0, tt] and m 2 + mj + m 2 3 = 1, (39) 

—%m\ m + im 3 y 

7.2 Extensions preserving parity 

The potential V(x), vanishing inside the box, is symmetric with respect to the point x = L/2. 
To make this symmetry explicit we shift the coordinate x to 

x 1 r 1 1. 

U — , «6 , -| — , 

L 2 L 2' 2 J ' 

and define 

V(u) = V(x), 4>e(u) = 4>e(x). 

In the new variable u the potential is even : V(—u) = V(u). It follows that, for a given energy, 
the eigenfunctions 4>e(u) and <Pe(— u) are solutions of the same differential equation and we 
can choose linear combinations of definite parity 4>e(u) ± 4>e(—u). 

As was already the case in the discussion of time reversal invariance, this argument is wrong 
since it overlooks the possibility for the boundary conditions (|3T| ) to break parity. Note that 
this point is often forgotten in Quantum Mechanics textbooks : there, one generally finds that, 
as soon as the potential is symmetric, the solution of the Schrodinger equation is of definite 
parity. It should be clear that the boundary conditions are essential. A good example to 
think about is the finite square well. The wave functions of its bound states are subject to the 

boundary condition J \(p(x)\ 2 dx < 00. As this condition is symmetric, the wave functions do 

have a definite parity. This is not the case for the diffusion eigenfunctions, for which one has 
an incoming and reflected wave for x — ► — 00, while for x — >• +00 one has only a transmitted 
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wave. In this second case the symmetry between x and — x is broken by the very conditions 
which characterize a diffusion experiment. 

We will therefore define parity preserving extensions of the Hamiltonian Hjj as the ones 
for which the eigenf unctions <Pe{u) verify 

\M-u)\ 2 = \Mu)\ 2 . (40) 

Here one finds (Appendix C) that all parity preserving extensions are given by m% = and so 
correspond to the matrix 

trriQ —m2 — imx \ 
tpe[0,n], m 2 + m 2 + m 2 = l. (41) 
m 2 — imi m / 



7.3 Extensions preserving positivity 

One of the most surprising facts, for a physicist, is the appearance of extensions with negative 
energies (these can be determined explicitly in some particular cases, see appendix B). 

From a theorem proved in |Tj| theorem 16, vol. 2, p. 44] one knows that only a finite 
number of negative energies can appear and that the sum of their multiplicities is at most 2. 
However, the determination of the U matrices with no negative eigenvalues, involves lengthy 
graphical discussions of equation (j35|) , which are fairly tedious. 

A partial answer to this problem is offered by an interesting theorem due to von Neumann 
(see |], p. 97]). It states that if A is densely defined and closed, then A^A is self-adjoint (and 
obviously positive). 

Let us apply this result to the operator (P = —iD, V (P)) defined in Subsection 5.3, whose 
adjoint was (P, "D max (0, L)). It follows that the operator 

(P 2 ,IMP 2 )), V^P 2 ) = {0 G D max (0,L) with 0(0) = 0(L) = 0}, 

will be self-adjoint. It does correspond to the extension with U = 1. 

If we take for operator (P, P max (0, L)), with adjoint (P, X> (P)), we are led to 

(P 2 ,P 2 (P 2 )), X? 2 (P 2 ) = {0 G V max (0,L) with 0'(O) = 0'(L) = 0}, 

a self-adjoint extension corresponding to U — —I. 

As a last example, we may start from (P,Ve), in which case von Neumann's theorem gives 
the self-adjoint extension 

(P 2 ,V 3 (P 2 )), V 3 (P 2 ) = {<peV m U0,L) with 0(L)=e^0(O), <//(L) = e»<//(0)}, 

corresponding to the matrix 




9 G [0,2tt]. 



As shown in the appendix B.2., for this choice of matrix U, the operators (P 2 , T>jj) and (P, T>g) 
have the same eigenf unctions. These extensions, (P 2 ,P[/) are really the square of the ones of 
the momentum operator (P, T>$) . 
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7.4 The infinite well as a limit of the finite one. 

Let us consider the standard problem of a particle of mass m in a one dimensional potential 
well of width L and depth Vq : 

V(x) = 0, xe)0,L[; V(x) = V >0, x &]0, L[. (42) 

A standard computation gives the bound states wave function 



, 2m(V - E) 



x < : <p n (x) = d n e px p , 

x>L : 4> n {x)= ±d n e-P^ (43) 



< x < L : <p n (x) — d n [cos kx + £ sin kx] 

with 

, k [2 



2m£ 



pVL v /(i + 2 /(pL))(l + F/p2) 

The positive integer n labels the (finite for a given value of Vo) family of solutions of the 
transcendental equation : 

and the ± corresponds to the (opposite) parity of the stationary state n, and to the relation 

cos(kL) + j sin(kL) = ±1 . 



When Vo is large, one finds for the spectrum (p ~ oo , Vq = y 2m ^ L2 ^> 1 ; k fixed) : 

k n L~Tm(l-2/v ) , E n ~E™{l-4/v ) (44) 
where the E 1 ^ 's are the infinite well energy levels (§), and for the stationary states : 

<p n (x < o) ~ ex P~ v o\ x / L \ ~° 

n (x > L) ~ ±\ t(— )exp-v (x/L-l) ~0 (45) 



L \v J 
<f) n {0< x < L) ~ v^TZ 



sinn7r h I — ) [cosn7r sinn7T— ] 

L \v J L nix L 



In that (fixed energy) infinite limit of the finite well, we see that the standard boundary 
conditions 0(0) = (f>(L) = are recovered. One could have considered a non-symmetric 
potential well such that V(x) = Vo for x < and V(x) = V\ for x > L with Vo ^ V\. 
Taking the limits Vo — > oo and V\ — > oo independently, leads to the same conclusions as for 
the symmetric case Vo = V\ considered here. 

This result is hardly a surprise since for fixed Vo we impose from the beginning the 
continuity of the wave function and its first derivative at x = and x = L. The wave func- 
tion in the classically forbidden region ( x < and x > L) is exponentially decreasing and is 
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damped off to zero in the Vo — ► oo limit. Combined with the continuity of <fi n {x) at the points 
x = and x = L this leads to 0(0) = (j){L) = (notice that in that limit the continuity of 
the first derivative of the wave function is lost). 

In many textbooks || vol. 1, p. 78], exercise 6.7, p. 396], this limiting process is argued 
to select the "right" boundary conditions for the self-adjoint extension of the Hamiltonian. 
In the same spirit, it would be tempting to consider the semi-axis case as a limit of a step 
potential. This selects uniquely the self-adjoint extension of the Hamiltonian such that 0(0) = 
(Subsection 6.2). However, for any finite height, the momentum P x has a unique self-adjoint 
extension, while for an infinite height, P x has no self-adjoint extension at all (see Subsection 
5.2)! 

This discussion shows that an infinite potential cannot be simply described by the limit of 
a finite one. 



8 Concluding remarks 

The aim of this article was twofold : first to popularize the theory of self-adjoint extensions 
of operators among people learning and (or) teaching quantum mechanics and second to point 
out some physical consequences which could be checked by experiment. 

For example the new spectra for a particle in a box should lead to different low temperature 
behaviours of the specific heat, following the lines of fT3"] , [0. Similarly, the boundary effects 
computed in || should be examined anew. 

Certainly the examples considered here are too simple, and are of questionable practical 
feasability. Our hope is that people will extend our analysis to the differential operators acting 
in three dimensional space which could lead to more realistic physical situations and put to light 
new phenomena : these developements could initiate the "physics of self-adjoint extensions" . 

Moreover, as previously seen, an infinite potential cannot be simply described by the limit 
of a finite one. This enforces interest in the large class of self-adjoint extensions described in 
this work : they deserve further study since they are all on an equal footing with respect to 
the principles of quantum mechanics. 

We have also emphasized in the previous Section the role of the symmetry properties (resp. 
reality properties) of the boundary conditions when the potential has some symmetry properties 
(resp. reality properties). Moreover, in subsection (5.3) we show that, in presence of an infinite 
discontinuity of the potential, the continuity of the wave function does not result from the 
principles of quantum mechanics. 

Last, but not least, let us mention other difficult problems which are not thoroughly dealt 
with in the standard teaching of quantum mechanics : the definition of higher powers of opera- 
tors ( to say nothing of their exponential !) and their commutators. This item was encountered 
in Section 2, where it was observed that H 2 is not the square of the operator H . On the 
contrary, in subsection 7.3, we have exhibited a specific extension of P 2 which is really the 
square of the extension Pq of P. 
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A Self-adjoint extensions of the momentum operator 

Let us consider the Hilbert space 7i = L 2 (a,b) . The maximal domain on which the operator 
P = —iHD has a well defined action has been called in Section. 5 T> max (a,b). It is the linear 
space of functions ip(x) constrained by : 

1. ip{x) is absolutely continuous |Tj| on [a, b]. 

2. ift(x) and ^'(x) belong to L 2 (a,b). 
It is useful to introduce the quantity 

Bty,® = \ % WA) - {^,P<f>)) = 5f0(6)0(6) -?(a)0(o)]- (46) 

A.l The operator P on the whole real axis 

The Hilbert space is Ji = L 2 (R) and the maximal domain of P is T> max (M). 
One can prove that for any i/j in this maximal domain, one has : 

lim ip(x) = 0. 

x— >±oo 

Note that this statement would not be true under the single hypothesis ip G L 2 (R). The 
symmetry of P is then, for 0, ip G £> max (IR), an obvious consequence of (f46"D- To prove that 
(P, Pmax(R)) is indeed self-adjoint, one should show that, if G L 2 (M) is such that 



G P max (M), 



ip'(x)4>(x) dx 



+oo \ 1/2 



then belongs to T , max (M). But it is easier to check this using von Neumann's theorem, which 
was done in Subsection 5.1. We have proven that the deficiency indices are (0, 0) and concluded 
that the operator (P, P max (K.)) is the unique self-adjoint extension of D. 

A. 2 The operator P on the positive semi-axis 

The Hilbert space is 7i = L 2 (0, +oo) and we take as domain 

X?o(P) = G P max (0, +oo) and ^(0) = 0}. (47) 
As in the previous subsection one can prove that lim ip(x) = 0. Then the symmetry of the 

x— >+oo 

operator P on T> (P) follows again from relation fllTf). 
The adjoint of (P,V (P)) is given by 

(pt = P, D(pt)=P max (0,+oo)). 

The double adjoint is simply 

Ptt = p, P(ptt)=p (p), 

which shows that (P,Vq(P)) is closed. 

However, as we checked in Subsection 5.2, the deficiency indices are (1, 0) and therefore, by 
von Neumann's theorem, (P, P(P)) has no self-adjoint extension. 
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A. 3 The operator P on a finite interval 

The Hilbert space is now Ti = L 2 (0, L) and we take 

P = -ihD, V (P) = {iJE V max (0, L), ^(0) = $(L) = 0}. 
The symmetry of P on T> (P) follows again from relation (f|Bp. Its adjoint is 

(p f = p, x?(pt) = p max ( ,D). 

Let us notice that the adjoint of (P, P(pt)) is (P,V (P)) which implies its closedness. 

In Subsection 5.3, we have obtained the deficiency indices (1, 1) and, from von Neumann's 
theorem, we know that the self-adjoint extensions are parametrized by U(l) i.e. a phase e tB . 

A. 4 Remarks 

1. In all cases we observe that the adjoint [P\ P(P')) has for domain P(P') = D max 
which is the largest domain in TC in which —ihD is defined. It follows that the actual compu- 
tation of the deficiency indices is always an easy task. 

2. Let us observe that for symmetric operators one has the hierarchy 

(P,X?(P)) C (P\V(P*)) 

which means that the adjoint is the "biggest". When self-adjoint extensions (P,T>q) do exist 
they must lie in the in-between, according to the scheme 

(P,X?o(P)) C (P,V e ) C {P\V{P^)) 

3. For further use, let us point out the useful theorem, proved in vol. 2, p. 90], stating 
that for a differential operator of order n with deficiency indices (n, n) all of its self-adjoint 
extensions have a discrete spectrum. 



B The spectra of the Hamiltonian in a box 

Starting from the boundary conditions ( |3~T| ) we now derive the equations giving the eigenvalues 
for all the extensions Hy. 

Let us consider the positive spectrum, the zero and negative ones being obtained in the 
same way and, as a matter of fact, obtained by substitutions as indicated in (|35|). 

s 2 

Denoting by E = — , with s > 0, the eigenvalues of Hu, and its eigenfunctions by 
L 2 



Ae isx/L + Be~ isx/L , $ = , (4 




one can easily check the relations 



L0'(O)-i0(O) \ / L0'(O)+i0(O) \ 

= i£(s)&, \=iM(s)$, (49) 

L6'(L)+i<b(L) / I L6'(L)-i6(L) / 
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with the matrices 

s 



Cs) 



M(s) 



1 -s - 1 

s + l)e is -(s-l)e- i! 
The determinants of these matrices are given by 
det M(s) = 2[i(s 2 + 1) sins - 2s cos s], 



3 + 1 

(s - l)e* 



-s + 1 

(s + l)e 



(50) 



det C(s) = -det M{s), 



from which it follows that C(s) and A4(s) have vanishing determinant if and only if s = 0. 
Using these notations the equations for the eigenf unctions become 

(C(s)-UM(s))$ = 0, (51) 

and for the spectra 

det (C(s) -UM(s)) = 0. (52) 

To get a more explicit form of the eigenvalue equation let us use some simple relations valid 
for arbitrary 2x2 matrices 

2 det A = (tiA) 2 -ti(A 2 ) det(A-B) = det A + det B + ti (AB) -trA-tvB. (53) 

For s we can write relation flo"2"| ) as 

det(£(s)7W~ 1 (s) - U) = 0, (54) 

where the matrix CM.^ 1 has the simple form 



C(s)M-\s) 



i\s 



1) sin s 
-2s lis 



-2s 
1 l)sins 



det M{s) 

Subsequent use of (|5^) in relation (p4]) and simple computations lead to 

2s [cos s(l - det £7) - tr (C/n)] + i sin s[(s 2 + 1)(1 + det U) - (s 2 - 1) tr U) = 0, 

valid for the positive non-zero spectrum. 

The parametrization of the matrix U given by 



(55) 



2s 



sin ip cos s + — tr (Mti) 



sin s 



), fl34]) simplifies relation fl55|) to 

> 2 + 1) cos?/> - -(s 2 - l)tr (M) 
2 



(56) 



a writing which exhibits the reality of the eigenvalue equation. It also displays a nice invariance 
under the transformation 



M 



2/ 



6 e [0,271-], 



(57) 



as it leaves trM and tr(Mn) unchanged. Let us point out that this invariance is specific of 
the spectra, not of the eigenfunctions. 

The strictly positive spectrum is then given by 

„2 



2s [sinip cos s — mi] = sins [cosip(s 2 + 1) — m (s 2 — 1)] 



E = — . 

L 2 



(5* 



The invariance (0) explains why the spectrum does not depend either of m 2 or of m 3 . 

An explicit solution of the eigenvalue equation (|58f) is clearly hopeless for the most general 
unitary matrix U. Nevertheless there are many special cases for which this can be achieved 
explicitly. We therefore classify the spectra as : 
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1. "Simple" if the eigenvalue equation can be solved explicitly. This happens for two 
families : 

mi = sinip = 0, or mo = costp = 0. 

2. "Generic" if this is not the case. Typically the "generic" spectra are solutions of at least 
one transcendental equation and only their large n behaviour can be obtained explicitly. 

B.l First family of "simple" spectra 

This first family corresponds to ip = and mi = and its matrix U has the form 

(m — im 3 —m 2 \ 
with m\ + m\ + m\ = \ <^> meS 2 . (59) 
m 2 m + im 3 J 

The eigenvalue equation reduces to 

sin s [(1 — m )s 2 + 1 + m ] = 0. 

Since m G [—1, +1] the factor in front of the sine never vanishes, so we get for spectrum 

s n = nn, n = 1, 2, . . . 

Note that the zero spectrum is easily checked to appear only for the extension with U — —I, 
while the strictly negative spectrum is given by 

sinhr [(m — l)r 2 + m + l] =0, 

which has always a solution, except for mo = ±1. We conclude to the negative energy 

2 1 + TP 1 1 + m , r 

r =- , — ► E = -— , m Q E -1,+1. 

1 — mo I/ 1 — mo 

Remark : In this family two and only two extensions (with mo = ±1) are therefore distin- 
guished by the absence of negative energies in their spectra. The first one is 



U = I 



s„ = nir, n = 1, 2 



0(0) = 4>{L) = I n ( x ) = sin . 

This is the "standard" self-adjoint extension considered in the textbooks on quantum mechanics 
§ p. 109], © p. 300]. 
The second one is 



U — —I 



717T, n — 0, 1, 



n'(0) = nU) = | n (a;) = 4/^cos(7i7r|) 

A different understanding of the absence of negative energies for these two extensions is 
given, using von Neumann theorem, in Subsection 7.3. 
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B.2 Second family of "simple" spectra. 

This second family corresponds to cosip = 0, or equivalently ip = tt/2, and m = 0. The 
corresponding matrix U is 

(m 3 m% — im 2 \ 
, with mj + ml + mj = 1. (61) 
mi + im 2 —m?, J 

Relation fl58|) reduces to 

cos s = m\. 

From (|6~TD we know that m\ G [—1, +1]. Excluding the values m\ = ±1, discussed in the final 
remark, the positive spectrum is 

+ cos -1 (mi) + 2nii, n = 0, 1, . . . 

cos" 1 (1) = tt/2. 

— cos (mi) + 2nir, n = l,2, ... 

As already observed, these eigenvalues are independent of m 2 and m 3 , but this degeneracy 
affects only the spectra, not the eigenfunctions. 
Let us observe that for the particular case 

/ e- i0 \ 
V e l6 J 

we have the full spectrum and eigenfunctions 

s n = 6 + 2nvr, <f> n {x) = - e 2 Mn+e/2*)x/L ? n = l,2,... 

L 

s n = -9 + 2nn, (j) n (x) = e ^i(-n+e/2^x/L^ n = 0,l,... 

Lj 

The exceptional cases 8 = and 9 = it are discussed in the next remark. The important point 
is that these eigenfunctions of P 2 are the same as for (P,T>g) given in Section 5.4 

Note that the zero spectrum is easily checked to appear only for the extension with U = T\ , 
while the strictly negative spectrum given by 

cosh r = mi 

is absent because from ( |B"TD we know that m\ e [—1, +1]. 

Remark : two extensions are distinguished by their doubly degenerate spectra. The first one 
corresponds to the periodic boundary conditions (the degeneracy of the energy s n is denoted 
by g n .) 

U = n ( s n = 2mr, n = 1, 2, . . . g n = 2, 

0(0) = 0(L), 0'(O) = d)'(L) ~^ { s = 0, 0o = 1, 
and the second one to the antiperiodic boundary conditions 

U = -n 

— ► s n = (2n + 1)tt, n= 0,1,2,... ^ n = 2. 

0(0) = -0(L), 0'(O) = -0'(L) 

The periodic boundary conditions may have a physical interpretation for rotational degrees of 
freedom of molecules H, vol. 2, p. 1202]. 
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B.3 The "generic" spectra 

We now exclude from our analysis the extensions with "simple" spectra. Switching to the 
variable t = tan -, the eigenvalue equation (|58j) becomes 

s [ , t \ 

'mi + sm^)t 2 H — [cosip{s 2 + 1) — m (s 2 — 1)] + rti\ — sin^ > = 0. (62) 



t 2 + l 

The overall factor should not be overlooked since it may vanish for t = oo. 

1 + t 2 3 
We organise the discussion of the "generic" spectra by distinguishing three different cases 

• mi = — sin ip 7^ 0. 

In this case the spectrum is 



cot - = — > s = (2n + l)ir, n = 0,1,2,... 

s (m — cos ib)s 2 — (m + cos ib) 

cot - = , s > 0. 

2 2s smip 



(63) 



Notice that smip cannot vanish (because then mi = and we are back to the first family of 
"simple" spectra). The numerator vanishes identically only for the second family of "simple" 
spectra, so we conclude that equation ( |63"D gives only "generic" spectra. 

• mi = sin ijj 7^ 0, 
in which case we have 

s 

tan - = — > s = 2nir, n = 1, 2, . . . , 
2 ' ' ' 

s (m — cos tb)s 2 — (m + cos tb) _ 

tan - = ■ , s > 0. 

2 2s smyj 

By the same argument as before neither the numerator nor the denominator can vanish, there- 
fore equation ( |5l| ) does give "generic" spectra. 

• mi ± smip 7^ 0, 

in which case the discriminant of equation (|62"D can be written 



A(s) = [(mo — cosip)s 2 — (m + cos^)] 2 + As 2 {m\ + ml), 

and is strictly positive because of the first term squared (otherwise we are back to the second 
family of "simple" spectra). 
The roots of 

tan £ = 1 Um a - m *)* - (m + cos,,) ± , (65) 

2 2s(m\ + smip) I J 

give spectra which are certainly "generic" . 

The equations giving the zero and the strictly negative spectrum can also be deduced as 
already explained in Section 6. 
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Let us conclude with a simple choice for the eigenfunctions 

A(s) = a(s-l) + [ 7 e- ls -l](s + l), fa 7 \ 

17 = . (66) 

7?(s) = a(s + l) + [ 7 e ls -l](s-l) = -A(-s), \ I 3 6 J 

This gives the relations 

<f){s; x) = A(s) e isx/L - A{-s) e ~ isx/L , 0(-s; x) = -<f){s; x). (67) 
The simultaneous vanishing of A and B signals a doubly degenerate spectrum. 

C Extensions preserving parity 

The eigenfunctions (^) write in u variable : 

j> E (u) = <j> E {x) = A{s)e ls/2 e isu + B{s)e~ is/2 e~ isu . 
Imposing the constraint (f40D gives 

lm(A(s)B(s)e is ) = 0. (68) 

It is important to observe that this relation should hold only when we take for s the actual 
spectrum given by relation (|58|) . 

Using for A(s) and B(s) the expressions given by (|66|), and after some algebra, one reduces 
the constraint (|68f) to 



(69) 



2s [(sin ip mo — cos ip m^) cos s — m^mi — m2m^\ = 

sins [(cos ip mo + sinip m^) (s 2 + 1) — (itlq + m 2 )(s 2 — 1)] . 

The m dependent terms disappear, thanks to relation (|58"D , and we are left with 

m 3 sin-^ [2s [cos?/' coss + m 2 ] + sins [sin^s 2 + 1) — m 3 (s 2 — 1)] } = 0. (70) 

One can check, by enumeration of all the cases, that the coefficient between braces never 
vanishes for m 3 ^ 0. 

We conclude that all the parity preserving extensions are given by m 3 = . Q.E.D 
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